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ABSTRACT 

Nighttime sky transparency statistics on Mauna Kea are reported based on 
data from the Canada-France-Hawaii Telescope SkyProbe monitor. We focus 
on the period beginning with the start of MegaCam wide-field optical imager 
operations in 2003, and continuing for almost three years. Skies were clear enough 
to observe on 76% of those nights; attenuations were less than 0.2 magnitudes 
up to 60% of the time. An empirical model of cloud attenuation and duration 
is presented allowing us to further characterize the photometric conditions. This 
is a good fit to the SkyProbe data, and indicates that Mauna Kea skies are 
truly photometric (without cloud) an average of 56% of the time, with moderate 
seasonal variation. Continuous monitoring of transparency during the night is 
necessary to overcome fluctuations in attenuation due to thin cloud. 



1. Introduction 

Sky transparency is a fundamental parameter governing the operation of any optical 
observatory. Best telescope performance comes when skies are photometric - cloud-free - 
and in general, efficiency at a given site depends directly on the fraction of time that skies 
are clear enough to observe. As more telescopes employ adaptive optics (AO) facilities 
employing lasers, even thin cloud can have a significant impact on telescope operations. 
Attenuation of the beam from cirrus will not only reduce laser return power and affect AO 
wavefront-sensor performance, but refiected Rayleigh-backscatter light may be of concern to 
other nearby optical observatories. 

Several sources of sky clarity measurements are available for Mauna Kea. The longest 
historical records date back to the 1970s, consisting of observers' nightly notes on sky condi- 
tions in the logbooks of the University of Hawaii 88-inch, United Kingdom Infrared Telescope 
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(UKIRT), and the Canada- France-H awaii Telescope (CFHT). Quantitative measure ments 



have been made in "campaign mode" (iMcCord and Clarklll979l : iKrisciunas et al.lll987l ). But 



only in the last few years has there existed dedicated instruments for measuring sky clar- 
ity. Two of these ar e the Continuous Camera (CONCAM) operated by the Night Sky Live 



network since 2005 (jShamir and Nemirofj|2005l ) and the Thirty Meter Telescope (TMT) All 



Sky Camera (ASCA) installed in 2006. Another is the CFHT SkyProbe. The CONCAM 
and TMT ASCA are all-sky imagers, and as such are best used to provide a panoramic 
view of cloud cover above the summit. SkyProbe is specifically designed t o measure sky 



transparency as an input to telescope data reduction (ICuillandre et al.ll2002l ) 



SkyProbe is fixed to the side of the CFHT. While the dome is open, once every minute 
it images bright stars in a 35 square-degree field overlapping the telescope field of view, 
performs automated photometry, and returns a measurement of attenuation corrected for 
zenith angle. It derives an absolute zeropoint based on the Hipparcos' Tycho catalog, which 
provides a solid photometric reference; uncertainty is on the order of 0.01 mag. An excellent 
feature of this design is that this measurement is in the direction of the telescope, and as such 
can be directly compared to a concurrent zeropoint measurement of another instrument, for 
example the MegaCam imager. The consistency of the SkyProbe data allows the ensemble 
of measurements taken during a night, or over more extended time periods, to be used 
to estimate general properties of the atmosphere. Nearby weather measurements from the 
CFHT meteorological tower can be used to check for consistency. SkyProbe does not take 
data at all times though, as it is shut off when the dome is closed. 

One utility of SkyProbe, helpful for the observer in pruning data of poor quality, is 
the detection of clouds. Thick clouds, say 0.5 magnitudes in extinction or more, are easily 
recognized by sharp fluctuations in attenuation. Less obvious are thin clouds, especially as 
their attenuations approach that of the atmosphere itself. To elucidate this, it is useful to 
separate attenuation into air and cloud components as 

■A = v4air + ^cloud = ^air,l-Z^ + ^cloud,!^^; (1) 

where cloud thickness T is analogous to airmass Z. Like Z, it parameterizes the bulk 
quantity of absorber without making any assumptions about its distribution along the beam. 
Constants Aair,i and Adoud.i indicate attenuations for Z = 1 and T = 1. Separated in this 
way, it can be seen that if conditions are usually photometric, 

^air,l = A, (2) 

where A indicates the mode of A, and 

^ioud,i =A-A, (3) 
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where A is the mean. So clouds are detectable if their attenuations are greater than A or 
when T varies on timescales shorter than Z - or both. 

Under photometric skies, astronomers often fix their photometry to a zeropoint assuming 
an airmass Z = 1 aX zenith for the observatory, although that is not strictly true. The density 
of air fluctuates with barometric pressure and temperature, and the co ncentrations of w ater 



vapour and other atmospheric constituents such as aerosols (see, e.g. IChambersI l2005l . and 
references therein). To first order, however, airmass will increase linearly with barometric 
pressure p and relative humdity p, 

Z=pp/{pp). (4) 

A complication is that p measured on the ground does not necessarily correlate with the 
bulk of the atmosphere. One extreme case is fog. But barometric pressure is generally 
stable during a night. So a plausible empirical model - requiring input only from SkyProbe 
- is that A, calculated nightly, tracks the observed atmospheric attenuation. And an even 
simpler model avoids the details of fluctuations in airmass entirely. This assumes that highs in 
barometric pressure (or any other atmospheric parameter) will be balanced by corresponding 
lows, giving a Gaussian distribution with a long term average of 

' ^°°expf-i^^)dZ = l, (5) 



V27rcr Jo ^ 2(T 

where a is the standard deviation in nightly A, expressed as a percent. 

When clouds are present, they may go unnoticed if they are thin, persistent and of 
uniform thickness. One possible means of detecting thin cloud is motivated by the fractal 
properties of cirrus (and stratus) clouds, in that their density and lifetime are cor related 



over many orders of magnitude, governed by a tight power law (llvanoval 120001 . |2002| ). If so, 
the number of SkyProbe samples of attenuation Adoud - the distribution of attenuations - 
should follow 

A^(^cioud) oc exp(-Acioud/aAioud,i), (6) 

where a is a constant and the thickness of clouds is given by T = v4cioudMcioud,i- For 
uniformly timed samples, the number of samples of a given attenuation should then be 
directly proportional to the duration of a cloud, S, and hence 

5 = (5oexp(-T/a), (7) 

where = /cieari^obs is the duration of clear skies during an observation lasting tobs, and 
/clear IS the fractiou of time that is clear. This implicitly assumes that samples are taken 
whatever the cloudiness, that is, without bias towards pointing the telescope at clearer parts 
of the sky, for example. That seems reasonable for conditions so good that clouds are hard 
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to detect. We will simply assume for now that samples are unbiased, and see if that still 
holds as clouds become thicker. Because the average thickness of a cloud is defined to be 

f = alog(5o/5) = l, (8) 

and since 

6/6o = 7^ exp ( - -) dT ^ exp(-l/«) (9) 

J max Jo ^ " ^ 

for maximal thickness Tmax > the thickness of a cloud would be given by 

T = l-alog(5/^)- (10) 

This suggests that an observer, after subtracting attenuation due to air, and knowing only 
a, could at least estimate the possible variation in cloud thickness and attenuation during a 
night, even if individual clouds themselves are hard to detect. 

In Section [2] we present the SkyProbe data and compare them to contemporaneous 
MegaCam photometric zeropoint measurements. Then, in Section [3] we show that our sim- 
ple models of Z and T are a good match to these data, reproducing both the effects of 
atmospheric attenuation and up to 3 magnitudes of extinction by clouds. A summary of the 
results follows in Section HI 



2. CFHT Data 

SkyProbe V^-band data were obtained from the period 1 January 2003 through 31 Oc- 
tober 2005. Each record runs from 6 PM to 6 AM Hawaii Standard time, with very few 
nights lost to the instrument being down. The first and last hour of data were deleted to 
avoid high sky brightness, which can lead to incorrect attenuation measurements due to 
slight nonlinearity in the SkyProbe detector. Weather data from the CFHT meteorological 
tower were obtained for the same period. The fraction of time when either SkyProbe did not 
take data, or the dome shutter was closed (including due to clouds) was 31.2%. It should 
be noted that the years 2003 through 2005 include some of the worst winter weather on 
record at CFHT. Even so, the mean attenuation from SkyProbe is 0.160 ±0.005 magnitudes 
per airmass, which is comparable, and even slightly better than the value of 0.206 ± 0.023 



reported by iKrisciunas et al.l (119871 ) based on UKIRT data spanning 1980 to 1986. This re- 
sult is also in good agreement with observer estimates of the sky-transparency for MegaCam 
during the period of overlap with the SkyProbe data. Tabled] is a breakdown of the average 
nights per year lost to each instrument due to weather and telescope downtime, along with 
the remainder with attenuation less than 0.2 magnitudes. 
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3. Analysis and Discussion 



Nightly modal SkyProbe attenuations and MegaCam zeropoints are plotted in Figured! 
along with air pressure, temperature, and relative humidity for the CFHT weather tower for 
the same period. Note the strong seasonal variation in air pressure and temperature, with 
winters clearly evident. There is no obvious seasonal dependence for relative humidity, apart 
from increased occurrences of saturation in winter. Modal attenuations from SkyProbe vary 
slowly from night to night, for the most part, with smooth oscillations over timescales of a 
month or more. The value of Aair,i inferred is 0.065 ± 0.027, taken to be the global mode of 
attenuations and the standard deviation of nightly modes over the same period. 

For comparison, equation 4 is overplotted on the SkyProbe data as a thick curve, boxcar 
smoothed to better show trends over a month. This agrees reasonably well with SkyProbe 
apart from some deviations, notably an underestimate for about three months in summer 
2005 (beginning with night 800). However, there are several sources of atmospheric atten- 
uation unaccounted for in equation 4. For example, Rayleigh scattering and absorption by 
ozone have a minimum of about 0.03 magnitudes at 0. 55 /xm for the altitude an d latitude 
of Mauna Kea, according to the standard calculation of iHayes and LathamI (119751 ). Another 
significant source of atmospheric attenuation is aerosols. Nephelometer measurements of 
aerosol absorption at nearby Mauna Loa Observatory were downloaded from the National 
Oceanic and Atmospheric Administration public archives. Averaged monthly absorptions 
(in the 0.5 fim channel) are plotted as a thin curve in Figured] for the same period. These 
have been normalized to a mean attenuation of 0.065 magnitudes, to help show trends rel- 
ative to equation 4. Variation in aerosol concentration is seen, and one period of elevated 
levels does correspond to the "enhancement" in SkyProbe attenuations in spring/summer 
2005 already mentioned. Another rise during the spring of 2003 (beginning approximately 
on night 50) might be evident in SkyProbe data; a third possibly in spring 2004 (night 450). 
Seasonal variation in th e concentration of aerosols has been seen at other observato ries, no- 
tably in central Europe (IReimann. Ossenkopf. fc Beyersdorferlll992l : iPakstiend l200ll ) . There 
the variation is thought to be the effect of higher summer temperatures raising smoke and 
soot into the upper atmosphere, combined with greater absolute humidity. It is not clear 
that the same mechanisms would be at work on Mauna Kea, although the seasonal pattern 
seems similar. One final "contaminant" of SkyProbe atmospheric attenuation estimates is 
evidently thin cloud itself. A few nights even have modal attenuations of up to 0.45 magni- 
tudes. Further color information might be helpful in discriminating between thin cloud and 
particulates, and a blue filter (to ma tch the Hipparchos catalog) has recently been installed 
on SkyProbe JCuillandre et al.lboosh . 



The bottom panel of Figure dl shows MegaCam zeropoints (subtracted from the mini- 
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mum zeropoint to indicate attenuation) over the same period. Overplotted are the nightly 
mode (for A < 0.2 mag) and median of SkyProbe, the latter highlighting cloudy periods, 
smoothed to help show trends over a month. Note that MegaCam zeropoints are already 
corrected by the observer for atmospheric attenuation, based partly on comparison with 
SkyProbe. It is possible that at times thin cloud may make this task more difficult, and the 
result more noisy, although the "striping" of zeropoints is due to a known processing problem 
unrelated to SkyProbe: residual uncorrected chip-to-chip offsets in MegaCam. The value of 
A inferred from MegaCam is 0.074 ± 0.031, in agreement with SkyProbe (0.065 ± 0.027). 

3.1. Separating Cloud from Air 

Figure [2] is a histogram of the SkyProbe data, plotted as attenuation in magnitudes. 
There is a strong peak at 0.065 magnitudes, indicating the global mode of atmospheric 
attenuation. To help see this, the distribution of nightly modal attenuations is overplotted 
as a dashed line. The second "bump" near 0.6 magnitudes is an instrumental artifact. It 
is caused by untracked exposures at the beginning of the night. At this time the telescope 
is sometimes parked at zenith and the resulting star trails are not processed properly by 
the automated software. A total of 145 nights identified with the resulting bogus zeropoint 
have been isolated and are not included in the analysis, although it is possible that some 
untracked exposures remain. 

Away from the peak associated with the zeropoint, a power-law behavior dominates. A 
linear least-squares fit to data with between 0.4 and 1 magnitude of extinction (solid line; 
extended above and below the fitted region) gives a power-law slope of -1.84. Comparing 
this with equation 6 gives a = — l/(slope x Ai^cioud) = — 1/(— 1.84 x 0.10) = 5.43. Note 
how well this curve fits - with constant a - even for attenuations approaching 3 magnitudes. 
We are not interested in this end of the distribution of course, as there is so little power 
there, but it is still worth mentioning. A bias against observing high A values - closing the 
dome for clouds, or systematically pointing the telescope to clearer parts of the sky - would 
presumably steepen the slope of this curve, perhaps even resulting in a knee at an attenuation 
where clouds become obvious, say, at 0.5 mag. A bend in the opposite sense would indicate 
observations taken preferentially in cloudy conditions. A more plausible scenario, consistent 
with the CFHT nightly log is that the dome is generally opened for nights that begin as 
mostly clear (with perhaps some cirrus), and observations continue as long as interruptions 
by thick cloud are brief. The dome is not opened during poor nights, as this leads to 
inefficient observing, and instead these nights are allocated to engineering tasks, if possible. 

The fraction of power under the curve (lightly shaded region) compared to that in the 
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total distribution is 18.6%. This is the total fraction of dome-open time that was cloudy. 
This indicates that for 81.4% of the SkyProbe samples, attenuation was due solely to air. 
SkyProbe was essentially always operational during the period discussed here and took data 
whenever the dome shutter was open. Taking this to be 68.8% of the time, conditions were 
photometric 0.814 x 0.688 = 56.0% of the time. The true value must be something more 
than this because downtime due to non-weather reasons (when the sky may not be cloudy) 
have been included, lumped in with times when the dome was closed (darkly shaded region). 
For example, because SkyProbe is attached to the telescope, datataking is subject to the 
telescope slewing to a new target. 

A more detailed histogram is shown in Figure [3l now plotted only for attenuations 
less than 0.8 magnitudes; the middle panel is that same, except with the nightly modal 
attenuation subtracted; and the bottom panel is a cumulative histogram of separated air 
and cloud. The thick black curve is the result of combining the distribution of nightly 
modal attenuations (thick dashed curve) with that for the cloud model, indicated by the 
shaded region. This provides a good fit, although it seems that the tail of the nightly modal 
distribution includes some nights when thin cloud was prevalent. To see this, the cloud 
model combined with the simple Gaussian model of atmospheric attenuation is indicated by 
the thin black curve. The dotted line indicates Gaussian noise of 0.01 magnitudes about a 
fixed atmospheric attenuation of 0.065 mag, that is, cr = and no cloud. It is the addition 
of this noise which sufficiently broadens the distribution of atmospheric attenuation that it 
leads to the nonphysical result of attenuations approaching zero, even becoming negative. 
This effectively restricts the accuracy of subtraction of the modal attenuation - only a spike 
could be perfectly subtracted - which is the cause of the scatter in the cloud measurements 
below about 0.2 mag, and ultimately makes the estimate of the clear fraction uncertain to 
within perhaps 2%. 

To see if any more information can be gleaned from the distribution of SkyProbe at- 
tenuations, for example, if the clear fraction of time varies seasonally, this calcuation was 
repeated after binning the data by calendar month. The results are presented in Figure |H 
The solid line is the calculation of cloud-free fraction of time, with the lightly-shaded region 
indicating cloud; dark shading indicates the dome is closed. Note that CFHT is very efficient 
in maximimizing the amount of time that the dome is open when skies are clear; this fraction 
of time varies moderately over the year, with perhaps a peak in late summer approaching 
60%. This and the values of other parameters appearing in equations 1 through 6 are given 
in Table [1 
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3.2. Simulated SkyProbe Data 

So far, we have only fit a power law to the SkyProbe data as per equation 6 without need 
for an explicit assumption about the duration of cloud-free skies other than to bin the data 
by month. It is interesting to consider for what values of 5o equation 10 is a good fit to the 
data. A computer simulation was developed to help explore this. Seven artificial SkyProbe 
datasets, assuming 6o = 1, 2, 6, 12, 24, 72 and 168 hours, were generated as follows: An 
initial cloud duration was selected at random from between zero and 6o- If a duration of 
zero was selected, then a period of 6q had no attenuation. Otherwise, T was calculated via 
equation 10. A power-law slope of a = 5.43 was used. In one-minute timesteps until the end 
of this period all datapoints each had this value. Then a new random duration was selected 
and the process repeated for a number of samples similar to the observations. Poisson noise 
of 0.01 magnitudes was added to simulate instrumental uncertainty. Data lost due to the 
dome being shut and to overheads was simulated by masking out 31.2% of nights with the 
highest median attenuation, just as CFHT closes for cloudy nights. As can be expected, since 
the mean of T has been defined to be unity, each of these simulations has the same mean 
attenuation, and correspondingly the same mean duration of clouds, that is Adond = 0.10 
magnitudes and 5 = 2.14 hours. 

The results are shown in Figure along with the SkyProbe data, displayed as grey- 
scale images. The nightly modal attenuation has been subtracted from the data to provide 
an estimate of attenuation due only to cloud. White indicates high attenuation; the same 
stretch is used throughout. Each pixel represents the average of three SkyProbe samples. 
Each row is one night of data, with a full year along a column. Instrumental noise is included, 
so clear skies appear as a uniform grey; black indicates the dome is shut. Small values of 
6o, less than 24 hours, seem to reproduce the short, thick clouds seen in the data. Larger 
values of 6o provide longer "clouded out" periods, some lasting several days, also evident in 
the data. But larger values of 6q produce too few short, thick clouds. 

It may be that a single value of 6q is insufficient to reproduce both long-term periods 
of cloud and the appropriate variation in attenuation during a night. This can be seen in 
Figure [6l which is a plot of the standard deviation of attenuation as a function of various 
time intervals. This has been calculated by dividing up each night's SkyProbe data into 
blocks of a given interval, and calculating the standard deviation over each. This analysis is 
restricted to attenuations of 0.5 magnitudes or less, so conditions excluding obvious clouds. 
Shown are values for 1, 2, 5, 10, 15, 30, and 60 minutes, and then every hour up to 10 hours. 
Note the strong power-law behaviour, as expected from the cloud model. Thin black curves 
indicate the results for choices of shown in Figure [5] plus that of (5o = 1 hour; the thick 
black curve is the average of these. Although the average of the models does not exactly 
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reproduce the data, it is a reasonably good fit, especially for observational durations longer 
than an hour. That the average of the models underestimates variation in attenuation on 
timescales less than an hour suggests that 6^ might include shorter timescales than this. 
And so it seems that a refined distribution of 60 - possibly spread over all timescales - could 
produce a better fit. But determining this distribution is beyond the scope of this paper. 
One path to it may be to investigate similarities with the Fried parameter ro associated with 
optical turbulence, where strong self-similar behaviour is also evident. 

A final result, evident in Figure El is that thin clouds set a threshold for photometric 
accuracy. Even if atmospheric attenuation is subtracted perfectly, if clouds less than 0.5 mag 
in attenuation are not accounted for, observations lasting an hour or more will typically incur 
a 0.03 magnitude absolute error in photometry. This "cloud-induced" error is comparable 
to the standard deviation in all MegaCam zeropoints (0.031 mag). As one might expect, 
the error component due to thin cloud could be reduced by confining observations to better 
sky conditions. It is not shown in Figure [6], but further restricting the SkyProbe analysis to 
attenuations less than 0.2 magnitudes reduces the slope of the cloud-induced error, yielding a 
0.01 mag standard deviation after 1 hour. Of course, the available observing time under these 
better conditions is reduced too, down to 60% based on Table [H only 56% of the time will 
have no clouds at all. This serves as a reminder to observers to obtain frequent photometric 
calibration; further validation of SkyProbe as an instrument dedicated to measuring sky 
transparency. 



4. Summary 

We have presented an analysis of CFHT SkyProbe data covering almost three calendar 
years beginning in 2003. The results are in agreement with MegaCam data during the same 
period. The SkyProbe data show that the V-band attenuation had a modal value of 0.065 
magnitudes corresponding to 56% of time with attenuation due to an atmosphere without 
clouds, and that 59% of the time V-band attenuation was 0.20 magnitudes or less (including 
both air and thin cloud). These estimates are conservative, as they include small telescope 
overheads. Some modest seasonal variation, ~ 4%, is seen in the fraction of time that 
skies are clear, but ultimately the accuracy in that estimate is perhaps ±2%, hmited by the 
photometric uncertainty of SkyProbe. As these results are for open dome shutter time on 
CFHT, and account for overheads, they can be used as inputs to designing AO with lasers, 
for example. 

A power-law model of cloud attenuation and duration has been presented which is a 
good fit to the SkyProbe data. It is consistent both with the mean attenuation due to 
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thin clouds, and the appearance of brief thick clouds at times causing the "clouding out" of 
nights. It is hoped that this simple model will lead to further investigation of the self-similar 
nature of clouds on Mauna Kea and elsewhere, which could possibly lead to improvements 
in photometric accuracy at the telescopes. 

We thank Derrick Salmon and Konstantinos Vogiatzis for their assistance in obtaining 
data from the CFHT weather tower and Chris Pritchet for his MegaCam zeropoints, and 
thoughtful comments. Wc would like to thank Glen Herriot and Jean-Pierre Veran for 
helpful conversations, and the Thirty Meter Telescope Site- Testing Advisory Group for useful 
comments on an earlier internal report on these results. This work includes archival data 
obtained from the National Oceanic and Atmospheric Administration, Climate Monitoring 
and Diagnostics Laboratory, Aerosols Group. 
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Table 1. Breakdown of SkyProbe and MegaCam Observing Statistics in Nights per Year^ 







Unobserved 






Attenuation (mag) 


Instrument 


Dome shut 


Overheads^ 


Total 


Observed 


> 0.2 


< 0.2 


SkyProbe 
MegaCam 


87 (23.8) 
87 (23.8) 


27 (7.4) 
37 (10.2) 


114 (31.2) 
124 (34.0) 


251 (68.8) 
241 (66.0) 


35 (9.6) 
20 (5.5) 


216 (59.2) 
221 (60.5) 



^Averages for a total sample of 1035 nights with SkyProbe and 541 nights with MegaCam. 
Fractional nights have been rounded down. Percent of total nights in the respective samples 
are shown in brackets. 

^Includes time lost during telescope slews, detector readout, technical faults, and engi- 
neering downtime. 
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Table 2. Atmospheric Parameters Derived from SkyProbe^ 



^air,l 


^cloud,l 


/clear 


(mag) 


a (mag) 


a (percent) 


0.065 ± 0.002 


0.027 ± 0.002 0.095 ± 0.003 


5.43 ±0.13 56.0 ±3.6 



^Averages after binning by calendar month. Quoted uncertainties are 
the standard deviations in each parameter. 
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Fig. 1. — Plots of nightly median barometric pressure, temperature, and relative humdity 
from the CFHT weather tower. Below are nightly modal V^-band attenuations from SkyProbe 
and gf-band zeropoints from MegaCam for the same nights. Atmospheric attenuation is 
overplotted for SkyProbe assuming equation 4 (thick curve), along with normalized monthly 
absorption due to aerosols (thin curve) estimated from nephelometer readings at Mauna 
Loa Observatory. Overplotted for MegaCam are nightly SkyProbe modes (thick curve) and 
medians (thin curve). 
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Fig. 2. — A histogram of the SkyProbe data. The peak at 0.065 magnitudes corresponds 
to the global mode of atmospheric extinction. To help show this, the distribution of nightly 
modal attenuations is overplotted as a dashed curve. A hnear least-squares fit to the data 
with 0.40 magnitudes or more gives a power-law slope of -1.84 (thick black line) corresponding 
to a value of a = 5.43. See text for details. 
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Fig. 3. — A histogram of the data (top panel; squares) that with nightly modal attenuation 
subtracted (middle panel; diamonds) and cumulative histogram (bottom panel) with a power- 
law model for clouds indicated by the shaded region. The two models for atmospheric 
attenuation are plotted as dashed curves, and the combination of these with the cloud model 
as solid curves. The dotted line indicates the minimum distribution of attenuations for Z = 1 
and photometric scatter of 0.01 mag. The fraction of time during the night without cloud 
is 56.0%. This is evidently not strongly dependent on the model of atmospheric attenuation 
used, as both provide good fits to the data. 
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Fig. 4. — A plot of cloud-free fraction of skies for the SkyProbe data calculated for each 
calendar month (thick curve); lightly shaded region is cloudy. This is corrected for the 
unobserved fraction of time during the night, indicated by the darkly shaded region. All 
unobserved samples are assumed to be cloudy. 
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Fig. 5. — Images of the SkyProbe data and simulations. Each pixel represents the average 
of three SkyProbe samples; white indicating high attenuation with the same stretch used 
throughout. Noise is included, so clear periods appear grey; times when the dome is closed 
are black. One night comprises a row. A year of data is shown. Simulations are shown for 
6o of 2, 6, 12, 24, 72, and 168 hours. None reproduces both the many short thick clouds seen 
in the data as well as longer "clouded-out" periods of several days. 
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Fig. 6. — A plot of standard deviations in SkyProbc attenuation due to cloud for observations 
lasting various time intervals, from 1 minute up to 10 hours. The nightly mode of attenuation 
has been subtracted for each night. The results of the simulations are shown as thin black 
lines, for 6o of 1, 2, 6, 12, 24, 72 and 168 hours; the average of those is indicated by a thick 
black curve. It would seem a model incorporating an appropriate distribution of Sq could 
reproduce the variation seen in the data. 



